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Abstract. Transverse hadron spectra from proton-proton, proton-nucleus and 
nucleus-nucleus collisions from 2 AGeV to 21.3 ATeV are investigated within two 
independent transport approaches (HSD and UrQMD). For central Au+Au (Pb+Pb) 
collisions at energies above ~ 5 AGeV, the measured AT* transverse mass spectra 
have a larger inverse slope parameter than expected from the default calculations. The 
additional pressure - as suggested by lattice QCD calculations at finite quark chemical 
potential / i q and temperature T - might be generated by strong interactions in the early 
pre-hadronic/partonic phase of central Au+Au (Pb+Pb) collisions. This is supported 
by a non-monotonic energy dependence of /{ pt ) in the present transport model. 
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Introduction 

Recent lattice QCD (1QCD) calculations at vanishing quark chemical potential and finite 
temperature indicate critical energy densities for the formation of a quark-gluon plasma 
(QGP) of ~ 0.6-1 GeV/fm 3 p. Such energy densities might already be achieved at 
Alternating Gradient Synchrotron (AGS) energies of ~ 10 AGeV for central Au+Au 
collisions 0 0 H- According to 1QCD calculations at finite quark chemical potential 
Tq a rapid increase of the thermodynamic pressure P with temperature above the 
critical temperature T c for a cross over (or phase transition) to the QGP is expected. 

Following the previous study [SJ we speculate that partonic degrees of freedom might 
be responsible for this effect already at ~ 5 A-GeV. Our arguments here are based on a 
comparison of the thermodynamic parameters T and /xg extracted from the transport 
models in the central overlap regime of Au+Au collisions jZ| with the experimental 
systematics on chemical freeze-out configurations jHJ in the T, /xg plane. The solid 
line in Fig. P characterises the universal chemical freeze-out line from Cleymans et 
al. 0 and the full dots with error bars denote the ’experimental’ chemical freeze-out 
parameters - determined from the thermal model fits to the experimental particle ratios 
P. The various smaller symbols (in vertical sequence) represent temperatures T and 
chemical potentials /xg extracted from UrQMD 1.3 transport calculations in central 
Au+Au (Pb+Pb) collisions at a/s = 200 AGeV, A[ a b = 160, 40 and 11 A-GeV jZj as a 
function of the reaction time in the center-of-mass (from top to bottom). 

During the non-equilibrium phase (open symbols) the transport calculations show 
much higher temperatures (or energy densities) than the ’experimental’ chemical freeze- 
out configurations at all bombarding energies (> 11 A-GeV). These numbers are also 
higher than the tri-critical endpoints and phase boundary extracted from lattice QCD 
calculations by Karsch et al. (S) (large open circle) and Fodor and Katz (5) (star with 
horizontal error bar). Though the QCD lattice calculations differ substantially in the 
value of /xg for the critical endpoint, the critical temperature T c is close to 160 MeV in 
both calculations, while the energy density is of the order of 1 GeV/fm 3 or even below. 
This diagram shows that at RHIC energies one encounters more likely a cross-over 
between the different phases when stepping down in temperature during the expansion 
phase of the hot fireball. 

Indeed, a hardening of the measured transverse mass (m t ) spectra in central Au+Au 
collisions relative to pp interactions mm from AGS energies on is observed. This 
increase of the inverse slope parameter T is commonly attributed to strong collective 
flow, which is absent in the respective pp or pA collisions. It has been proposed pj2| to 
interpret the high and approximately constant K ^ slopes above ~ 30 AGeV - the ’step’ 
- as an indication of the phase transition. 

In this contribution we explore whether the pressure needed to generate a large 
collective flow to explain the hard slopes of the K ± spectra with a ’plateau’ at SPS 
energies is produced in the present transport models by interactions of hadrons or 
whether additional partonic contributions in the equation of state might be needed 


r Pfp f n qiiptqp PrpQQ'iirp nnrl Rirnn npn pqq Dimnmira in PpI nfini qH r Hpnnn Tnn Ppnriinn q 3 



M-jj [MeV] 

Figure 1 . Schematic phase diagram in the T — /j,b plane. The solid line characterises 
the universal chemical freeze-out line from Cleymans et al. [R; and the full dots (with 
error bars) denote the ’experimental’ chemical freeze-out parameters from Ref. [Hj- 
The various symbols represent temperatures T and chemical potentials /is extracted 
from UrQMD 1.3 transport calculations in central Au+Au (Pb+Pb) collisions at 21.3 
A-TeV, 160, 40 and 11 A-GeV fj) ( see text). The large open circle and the star indicate 
the tri-critical endpoints and phase boundary from lattice QCD calculations by Karsch 
et al. 0 and Fodor and Katz respectively. The horizontal line with error bars is 
the phase boundary from [jj. 


to explain these effects (for further details the reader is referred to t Bj). To understand 
whether a failure of the present models indeed hints a QGP onset, we explore two 
distinct effects that might result in a substantial increase of the transverse pressure: I) 
initial state Cronin enhancement and II) heavy resonance formation. 

The Models 

In our studies we use two independent relativistic transport models that employ hadronic 
and string degrees of freedom: UrQMD [HI US! and HSD pa IT7j . They take into 
account the formation and multiple rescattering of hadrons and dynamically describe the 
generation of pressure in the hadronic expansion phase. This involves also interactions 
of leading pre-hadrons that contain a valence quark (antiquark) from a primary ’hard’ 
collision (cf. Refs. JXE, EEHD- Note that, in these models, only hadrons, valence 
quarks and valence diquarks and their interactions are treated explicitly. Gluonic 
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degrees of freedom are not treated explicitly, but are implicitly present in strings. This 
simplified treatment is generally accepted to describe proton-proton and proton-nucleus 
interactions. Here we test whether this description is still valid for the more complicated 
nucleus-nucleus collisions, where large energy densities can be reached over extended 
volumes. 


Transverse Dynamics in Small and Large Systems 

Let us start by ’’benchmarking” the model calculations with pA data. Fig. [21 shows the 
results for the inverse slope parameters T for various reactions - see figure caption for 
details. It can be seen that the models reproduce the transverse slope parameters of 
different particles produced in pA interactions with targets from Be to Pb reasonably 
well. 




Figure 2. Inverse slope parameters T for tv ± ,K + and K~ at midrapidity from pA 
reactions at 14.6 GeV/c (A= Be, Al, Cu, Au) - left part and at 450 GeV/c (A= Be, 
S, Pb) - right part, from HSD (open symbols) and UrQMD 2.0 (closed symbols). The 
full symbols in the left part correspond to the midrapidity data {(yiab) = 1-5,1.7,1.9) 
from the E802 Collaboration 2Z], in the right part to the NA44 data |28| at 2.4 < 
yiab < 3.5, pt < 0.84 GeV/c for K + ,K~ and at 2.4 < yi a b < 3.0, pr = 0.3 4- 1.2 
GeV/c (full diamonds) and 3.1 < yi a b < 4.0, pr < 0.64 GeV/c (full squares) for tt + . 

We continue with nucleus-nucleus collisions, where Fig. El summarises our results: 
the dependence of the inverse slope parameter T on is shown and compared to (partly 
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preliminary) experimental data from [TUI 12m 2T1I221I2U1 'TJj for central Au+Au (Pb+Pb) 
collisions (l.li.s.) and [22, 1UU1IUU1 for pp collisions (r.h.s.). The upper and lower solid lines 
(with open circles) on the l.h.s. in Fig. 01 correspond to results from HSD calculations, 
where the upper and lower limits are due to fitting the slope T, an uncertainty in the 
repulsive A'^-pion potential or the possible effect of string overlaps. The solid lines with 
stars correspond to HSD calculations with the Cronin effect. The dashed lines with 
open triangles represent slope parameters from UrQMD 1.3, the dot-dashed lines with 
open inverted triangles correspond to UrQMD 2.0 results, which are well within the 
limits obtained from the different HSD calculations without the Cronin enhancement. 
The dotted lines with crosses show the UrQMD 2.1 results that incorporate high mass 
resonance states up to m R < 3 GeV. 

The slope parameters from pp collisions (r.h.s. in Fig. 0 are seen to increase 
smoothly with energy both in the experiment (full symbols) and in the HSD calculations 
(full lines with open circles). The UrQMD 1.3 results are shown as open triangles 
connected by a solid line and systematically lower than the slopes from HSD at all 
energies. When including jet production and fragmentation via PYTHIA in UrQMD 
2.0 (dot-dashed lines with open inverted triangles) the results become similar to HSD 
above a/s = 10 GeV demonstrating the importance of jets in pp reactions at high energy. 

Coming back to the slope parameters of Ii^ mesons for central Au+Au/Pb+Pb 
collisions (l.h.s. of Fig. 0 we find that the Cronin initial state enhancement indeed 
improves the description of the data at RHIC energies, however, does not give any 
sizeable enhancement at AGS energies. Here UrQMD 2.1 (dashed lines with crosses) 
with the high mass resonance states performs better: Including high mass resonances one 
obtains more reasonable results for K + mesons, however, fails by 10 to 15% for pions as 
well as anti-Kaons. In this context it is interesting to note that the experimental results 
on C+C and Si+Si at 158 AGeV show small slopes [2S] and are therefore in agreement 
with the models (1). 

What is the origin of the rapid increase of the slopes with energy for central Au+Au 
collisions at AGS energies and the constant slope at low SPS energies (the ’step’ in the 
Kaon temperature), which is missed in presently employed transport approaches? 

Elliptic Flow 

To disentangle the effects of high mass hadron states from a possible phase transition 
scenario we suggest to study the energy excitation function of the elliptic flow of pions 
(or negatively charged hadrons). The ( p R ) and v 2 excitation functions are depicted in 
Fig. 0(left). One clearly observes a monotonic rise in the mean transverse momentum of 
the pions with increasing energy. However, the elliptic flow behaves non-monotonic and 
shows a distinct maximum around 30-40 GeV beam energy. This phenomenon can be 
pronounced more clearly by the scaled elliptic flow (v 2 /(pr)) as shown in Fig. 01 (right). 
It was pointed out HE EDI that at high energies the differential v 2 { j > t ) of charged hadrons 
is approximately proportional to p R , such that the averaged v 2 oc ( pr )• In fact, when 
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Figure 3. Comparison of the inverse slope parameters T for K + and K~ mesons from 
central Au+Au (Pb+Pb) collisions (l.h.s.) and pp reactions (r.h.s.) at midrapidity as 
a function of the invariant energy y/s from HSD (upper and lower solid lines with 
open circles), UrQMD 1.3 (dashed lines with open triangles), UrQMD 2.0 (dot-dashed 
lines with open inverted triangles), UrQMD 2.1 (dotted lines with crosses) with the 
data from Refs. QU1 EDI EH E2 E3 El for AA and C3 E3 EDI for pp collisions The 
upper and lower solid lines in the left diagrams result from different limits of the HSD 
calculations as discussed in the text while the solid lines with stars correspond to HSD 
calculations with the Cronin initial state enhancement. 


divided by the average transverse momentum, the scaled elliptic flow in that energy 
regime becomes nearly constant mm, even though the measured v 2 increases from 
SPS energy to RHIC energy. 

To emphasise deviations from the natural scaling v 2 oc ( p?) , we plot the excitation 
function of v 2 /(pt ) in Fig. El (tight). At lower energies there is a systematic increase 
of v 2 relative to (pt) in the model and in the data. However, above the SPS energy 
regime, one clearly observes that the data for the scaled elliptic flow is constant and 
independent of energy while the model yields first a sharp decrease of the v 2 /(pt) which 
then levels-off at roughly half the experimentally observed value. 

The present non-equilibrium study however, suggest that the initial increase of the 
scaled elliptic flow up to SPS energies, might be due to viscosity effects (decreasing 
mean-free-path) in the hadronic gas. At higher energies, the predicted elliptic flow 
breaks down in the model calculation because of the increasing dominance of string 
dynamics. The measured data however, supports a hydro dynamical behaviour of the 
matter in the early stage with very small mean free paths. It should be noted that a 
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similar minimum in the elliptic flow was first predicted by hydrodynamics. In contrast 
to the scenario discussed here, this decrease of v 2 was associated with the softening of 
the equation of state in the phase transition region iEZUSS]. 
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Figure 4. Left: Calculated excitation functions of mean transverse momentum 
and elliptic flow of pions at midrapidity. Right: Excitation function of v 2 /{pr} of 
Pions in mid-central collisions from top AGS to RHIC energy from UrQMD (squares). 
Negatively charged particle data (diamonds) for v 2 are taken from Ref. m and (pr) 
from Refs. ISSIIMIETIEHI- 

Thus, from the lack of initial pressure we conclude that the system (at least at 
RHIC energies) seems to spend a considerable amount of time in the QGP phase with 
an equation of state harder than the employed hadron/string gas equation of state. 
This argument is well in line with the studies on elliptic flow at RHIC energies, which is 
underestimated by ~ 30% at midrapidity in HSD JJfSj and by a factor of ~ 2 in UrQMD 
1-3 M- It is our opinion that strong pre-hadronic/partonic interactions might cure this 
problem. 


Conclusion 

In conclusion, we have found that the inverse slope parameters T for K ± mesons from 
the HSD and UrQMD 1.3 transport models are practically independent of system size 
from pp up to central Pb+Pb collisions and show only a slight increase with collision 
energy. The calculated transverse mass spectra are in reasonable agreement with the 







Transverse Pressure and Strangeness Dynamics in Relativistic Heavy Ion Reactions 8 

experimental results for pp reactions at all bombarding energies investigated as well as 
central collisions of light nuclei (C+C and Si+Si) (cf. Ref. |6]). The rapid increase of 
the inverse slope parameters of Kaons for central collisions of heavy nuclei (Au+Au or 
Pb+Pb) found experimentally in the AGS energy range, however, is not reproduced by 
both models in their default version (see Fig. ED- 

We have discussed scenarios to improve the description of the experimental 
data. However, no fully convincing results could be obtained for all observables and 
bombarding energies simultaneously. 

From comparison to lattice QCD calculations at finite temperature and baryon 
chemical potential /ib from Refs. [Hj and j9jj as well as the experimental systematics 
in the chemical freeze-out parameters (cf. Fig. ^), we infer that the missing pressure 
above 30 GeV beam energy might be generated in the early phase of the collision by 
non-perturbative partonic interactions. However, to fully clarify this issue will require 
a systematic quantitative comparison with hydrodynamic models from the lowest AGS 
energy to the highest RHIC energy. 
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